Lack of Self Averaging and Finite Size Scaling in Critical Disordered Systems 
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We simulated site dilute Ising models in d = 3 dimensions for several lattice sizes L. For each 
L singular thermodynamic quantities X were measured at criticality and their distributions P{X) 
were determined, for ensembles of several thousand random samples. For L —* oo the width of P(X) 
tends to a universal constant, i.e. there is no self averaging. The width of the distribution of the 
sample dependent pseudocritical temperatures T c (i,L) scales as 5T C (L) ~ L~~ and not as ~ L~2 . 
Finite size scaling holds; the sample dependence of Xi(T c ) enters predominantly through T c (i,L). 
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Phase transitions in systems with quenched disorder 
are of considerable theoretical and experimental inter- 
est and have been the subject of intensive investigations. 
Nonetheless, there are only few general exact results on 
such systems and some of these raised questions which 
were left unanswered for more than a decade. For ex- 
ample, the Harris criterion jl[ states that the condition 
for the stability of a pure fixed point against disorder is 
l/^pure < d/2. This criterion was derived assuming that 
fluctuations of the local critical temperature T c j of a re- 
gion of size £f scaled as 5T C ~ £ _d / 2 . Another result by 
Chayes et al 0] stated that 1/v < d/2 for any disordered 
system. This result is identical with the Harris criterion 
for disordered systems governed by a pure fixed point. 
For systems governed by a random fixed point it raised 
the question whether the relation ST C ~ £~ d / 2 holds in 
this case as well? 

Another important and relevant issue concerns the 
meaning of measurements done on disordered systems at 
or near their transition points. All measurements (exper- 
imental and numerical) are obtained for finite systems; 
moreover, each such system constitutes a particular re- 
alization of the quenched randomness. Hence it is natu- 
ral to ask to what extent are results obtained for a sin- 
gle system representative of the general class to which 
it belongs? The answer to this question hinges on the 
important issue of self-averaging. If a quantity is not 
self averaging, increasing its size L does not improve the 
statistics of its measurement (sample to sample fluctua- 
tions remain large) . Whereas it has been known that for 
(spin and regular) glasses |S| there is no self-averaging 
in the ordered phase, the discovery that there is no 
self averaging for random ferromagnets at their critical 
point came as somewhat of a surprise. The numerical 
work on which this claim was based has been backed up 
with a finite-size scaling ansatz Q which seemed to fit the 
data quite well. Subsequently the issue has been investi- 
gated by Aharony and Harris || (AH), who provided a 
theoretical understanding of the absence of self averag- 
ing in critical random ferromagnets. The central point 
of the AH argument was that a random fixed point is 



characterized by a distribution of non-zero width of some 
measurable quantities (such as pseudo-critical tempera- 
tures, susceptibilities etc) and hence if an ensemble of 
systems flows to a such a fixed point, its corresponding 
properties must also be distributed in a similar way. The 
general argument was supported by a Renormalization 
Group calculation in d = 4 — e dimensions. The AH work 
led to several predictions, some of which were in clear dis- 
agreement with the scaling theory [Q. Since neither the 
general arguments nor results obtained by e-expansion 
can be viewed as definitive, independent confirmation is 
desired. Furthermore, recently claims were made j(| to 
the effect that the manner in which the sample to sample 
fluctuations of the pseudoritical temperature scales with 
size governs one's ability to observe the "true" critical 
exponents of random systems, quantum and classical. 

In order to investigate these fluctuations, the extent 
to which our finite size scaling ansatz holds and to con- 
firm or disprove the AH results we carried out exten- 
sive simulations of the 3 — d Ising model with site dilu- 
tion. Since for the pure model the specific heat exponent 
a p = 0.11>0[Q, randomness is relevant and the critical 
behavior is governed by a random fixed point. Various 
aspects of the critical properties of this model have been 
measured very carefully [^J. In particular, a, the spe- 
cific heat exponent of the random model is negative [0, 
consistently with || . 

To state the main results of our study we first re- 
capitulate a few definitions. We consider an ensemble 
of systems of linear size L, denoting a particular re- 
alization of the randomness by i. For each such sys- 
tem we measure, at temperature T, various thermody- 
namic densities Xi(T,L) (such as the susceptibility per 
site X — X, magnetization m, etc). Denote by T c (i,L) 
the pseudocritical temperatures of each of these samples 
(obtained, for example, by locating the maximum of \, 
see Fig 1). The sample- averaged pseudocritical temper- 
ature, T C (L) = [T c (i,L)] approaches, for L — ► oo, the 
asymptotic limit T C (L) — > T c . Our main results can be 
summarized as follows: 

1. Lack of self-averaging: The values of Xi measured 
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at T c for many samples of size L, are distributed 
with mean [x(£)] and variance V X (L). The nor- 
malized square width 



Rx = Vx/ix? 



C as L 



i.e. goes to a constant, for large L, as predicted by 
AH, in contradiction to our prediction R ~ L a ' v . 
The distributions of various thermodynamic quan- 
tities do not become sharp in the thermodynamic 
limit and the system is not self-averaging. 

2. Universality: AH predicted that the constant C is 
universal. To test this, we studied different dilu- 
tions and two random ensembles; one with a fixed 
concentration of magnetic sites (grand canonical) 
and the other with a fixed number of spins, placed 
at random on the lattice (canonical). Grand canon- 
ical ensembles with different concentrations indeed 
have the same asymptotic width, but two different 
ensembles seem to have different values 

R (even when the concentration in one equals the 
fraction of spins of the other). 

3. The pseudocritical temperatures T c {i, L) are dis- 
tributed with a width 5T C (L) which scales as 



ST C (L) ~ L~ 



(1) 



This behavior was sugested by AH, reasoning that 
when combined with the finite size scaling theory 
of ID it gave rise to the lack of self averaging dis- 
cussed in item |l|. However this reasoning depended 
on the correctness of the finite size scaling theory. 
Here we have directly shown the validity of (Q) inde- 
pendently of any additional assumptions. We rule 
out the possibility that 5T C (L) ~ L~ d / 2 as was as- 
sumed by us Q and listed by others (|] as "the most 
likely scenario" . The latter behavior occurs only in 
systems governed by a pure fixed point where the 
Harris criterion is valid. Note that the critical con- 
centration in percolation JllJ] behaves as in ([!]) (see 
also a similar phenomenon in the random field Ising 
model ||). 



4. The finite size scaling form assumed in 
X l {T,L)^LPQ l (t l L L «) 



(2) 



where ti is a sample-dependent reduced tempera- 
ture, 



t i = [T-T c (i,L)]/T c 



(3) 



is indeed correct. Particularly, the sample depen- 
dence of Xi is predominantly due to the sample 
dependence of ti and not due to the sample depen- 
dence of the function Qt. 



5. However, the scaling function Qi{x) does depend 
on the realization i. 

6. It may be computationally advantageous to mea- 
sure various quantities, for instance, to find critical 
exponents, at the sample-dependent pseudocritical 
point T c (i,L) (where ti = 0), rather than at T c . 
This is so since, as a consequence of item [| the 
variance at T c (i,L) is much smaller. 

We present now numerical evidence for each of the state- 
ments made above. 

The site- dilute Ising model: Each site of a cubic lattice 
is either occupied by an Ising spin or empty. Here we 
report results obtained for a grand-canonical ensemble 
of samples in which the occupation of each site was de- 
termined independently with probability p — 0.8. At 
this probability Heuer found the fastest crossover to 
the asymptotic (random) critical behavior ||. We used 
the Wolff [jl3| single cluster algorithm jL4} with skewed 
periodic boundary conditions jl5) on lattices of sizes 
L = 4, 8, 16, 32, 64. For each of these sizes we simulated, 
respectively, 

n L = 10000,4000,32000,4000,1479 

different random samples. High-precision measure- 
ments of various critical properties yielded T c — 
3.49921(3), a/v = -0.066(9), (3/v = 0.505(2), jju = 
1.990(4), l/v = 1.467(5), in good agreement with g. 

Lack of self -averaging, distribution of the critical sus- 
ceptibility: Fig 1. presents curves of x{T) vs temperature 
for two samples of size L = 32. 
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FIG. 1. Susceptibility per spin of two systems of size 
L — 32, generated with site occupation probability p = 0.8, 
versus temperature. T c (i,L) denotes the pseudocritical tem- 
perature of sample i, identified by the maximum of the sus- 
ceptibility. T c is the value approached by the mean of these 
in the L — > oo limit. Xi(T c ) are the values of x as measured 
at T c for sample i. 

The values of Xi(T c ) arid the pseudocritical tempera- 
tures of the two samples are clearly identified. From mea- 
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suring Xi(T c ) for an ensemble of samples we constructed 
the histogram of Fig 2, for various sizes L. 
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FIG. 2. Distribution of the critical susceptibilities, mea- 
sured at T c for different samples and different sizes L = 16 
(thin dotted line), L = 32 (thick dashed) and L = 64 (thin 
solid). For each L the values of Xi(T c ,L) were normalized by 
the ensemble average of all samples of that size. The inset 
shows the second moment of the distribution as a function of 



L; evidently 



C. 



As evident from the histograms and from the inset, 
the widths of these distibutions approach a constant for 
increasing L. This is the result of AH, whereas our scal- 
ing theory [Q would predict R x ~ L a l u . The source of 
the discrepancy is an assumption we made regarding the 
distribution of the pseudocritical temperatures; had we 
assumed that (^) holds, our scaling ansatz would have 
also produced the correct result R x ~const. Therefore 
confirmation of this result may serve as an indication of 
the validity of ([l]) as well. Nevertheless, since this inter- 
pretation depends on the validity of further assumptions 
implicit in our scaling theory, it is important to verify 
independently and directly scaling and (Q). 

Universality: We studied two more random ensembles: 
another grand canonical one, but with occupation proba- 
bility p = 0.6 and a canonical one with a fixed number of 
0.6L d spins, randomly placed in each sample. We found 
that the limiting values of both R x and R m are the same 
for the two grand-canonical ensembles (with p = 0.8, 0.6) 
indicating universality, but for the canonical ensemble 
different values were obtained. These results will be pre- 
sented and explained elsewhere ||. 

Distribution of the pseudocritical temperatures: We es- 
timated, using the histogram reweighting method [fl6l| , 
the pseudocritical temperature T c (i, L) of every sample 
of our ensembles. This was done in an iterative way; the 
first guess for T c (i, L) was T c - from data collected at this 
point we calculated the susceptibility as a function of T; 
the temperature at which it had its maximum was our 



next estimate for T c (i, L) where more data were collected 
and so on. The procedure converged in less than ten it- 
erations within the resolution imposed by the statistical 
error in determining the temperature of the maximal \. 
The resulting histograms for several sizes are collapsed 
on Fig 3 using the value of yt = - = 1.467(5) which was 
determined independently |1C|] by simulations at T c . 
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FIG. 3. Distribution of scaled pseudocritical temperatures 
obtained for sizes L = 16 (thin dotted line), L — 32 (thick 
dashed) and L = 64 (thin solid), using yt = 1/u = 1.467 and 
T c = 3.4992. 

The means and widths of these distributions are 

T C (L) « T c - 0.532 • iT 
ST C (L) PS 2.13 • L- X l v (4) 

The width of the distribution exceeds by a factor of 
nearly 4 the shift of its average from T c . Hence the large 
fluctuations in x(T c ) are due to the large value of 5T C (L); 
when we perform measurements at the fixed temperature 
T c , some samples we deal with will be considerably above 
their pseudocritical temperature and some below! (see 
Fig 1) We have also made straightforward fits of the vari- 
ance ST C (L) 2 to the form 5T C (L) 2 - L~ 2p , and fits of the 
shift of T C (L) to T C (L)-T C - r A ; we found p = 1.449(8) 
and A = 1.42(4). Clearly p is within errors of the shift 
exponent A and its error bar is small enough to conclude 
that p d/2 = 1.5. Since one expects that X — yt and 
we have also determined independently y t = 1.467(5) our 
results strongly suggest that all these exponents are ac- 
tuary the same p = A = 1/v. 

Testing the scaling ansatz: Eq. (^) represents an im- 
plicit assumption; that the most important effect of the 
sample dependence of each realization of the randomness 
can be absorbed in the pseudocritical temperature of the 
particular sample. Our first task was to check the extent 
to which this holds, i.e. to what extent can one collapse 
data obtained at different sizes and temperatures for dif- 
ferent samples? Fig. 4 presents the magnetization for 
two system sizes. The data collapse quite well onto two 
branches of a function, one below and one above T c (i, L), 
lending support to the validity of the scaling ansatz. We 
succeeded to fit thousands of data points at three differ- 
ent sizes to two scaling functions (corresponding to the 
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two branches). This may be interpreted as supporting a 
strong scaling hypothesis, e.g. that for large L 



Xi(T,L) fts L p Q(iiL») 



(5) 



Here we modified (g) by dropping the sample-dependence 
of the scaling function, i.e. assuming that the entire sam- 
ple dependence can be absorbed in ij. On the other hand, 
as evident from Fig. 4, the data points do exhibit consid- 
erable scatter about the main trend. To investigate the 
extent of violation of such a strong scaling relation, we 
studied the distribution of the susceptibility maxima. 
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FIG. 4. Scatter plot of the scaled magnetization, mea- 
sured at T c for ensembles of samples of two different sizes, 
plotted versus the scaled temperature U = {T c — T c (i, L)}/T c . 
The lines are the scaling functions fitted separately for three 
different sizes, showing remarkable agreement. The scatter 
of the points from the lines is due partially to measurement 
error (thermal) and partly to genuine sample dependence of 
the scaling function. 

Sample dependence of the scaling function: If, in- 
stead of %i(T c ,L), we measure for each sample and 
size the maximal value of the susceptibility, xT ax — 
x[T c (i,L),L)], we accumulate all our data at the same 
value of the scaling variable tiL~ = 0. Hence if for large 
L the scaling function goes to a sample independent form, 
the distribution of xT ax should approach a S function. 
This distribution is presented in Fig 5: even though it is 
much narrower (by a factor of m 70) than the distribu- 
tion of Fig 2, as the inset shows, its width also goes to a 
constant. Hence a strong scaling hypothesis such as (||) 
can be viewed only as an approximation. 

Consequences for efficient simulations: In order to ac- 
quire data for large systems with high precision (say to 
estimate exponents from finite size scaling analysis) the 
commonly accepted procedure is to perform, for all sizes, 
simulations at one temperature, T c . The error inherent in 
doing this, due to sample to sample fluctuations, is much 
larger than that of taking measurements at T c (i, L). The 



latter procedure involves a different computational over- 
head, that of determining T c (i, L). Clearly, this may be 
advantageous if the width of the quantities measured at 
the pseudocritical temperature of each sample is much 
smaller than that of the data taken at T r . 
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FIG. 5. Distribution of the normalized maximal suscep- 
tibilities for samples of sizes L = 16 (thin dotted line) and 
L — 32 (thick dashed). The horizontal scale is the same as 
Fig 2, to emphasize how small is the width of this distribution, 
which also goes to a constant (see inset). 
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